Avoided level crossing spectroscopy with dressed matter waves 
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We devise a method for probing resonances of macroscopic matter waves in shaken optical lattices 
by monitoring their response to slow parameter changes, and show that such resonances can be 
disabled by particular choices of the driving amplitude. The theoretical analysis of this scheme 
reveals far-reaching analogies between dressed atoms and time-periodically forced matter waves. 
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Recently it has been demonstrated experimentally that 
a macroscopic matter wave of ultracold bosonic atoms 
confined in an optical lattice can be controlled in a 
systematic manner by strong, off-resonant time-periodic 
forcing: Under suitably selected conditions, "shaking" 
the lattice with kilohcrtz frequencies mainly effectuates a 
modification of the tunneling matrix element connecting 
adjacent lattice sites. In the regime of weak interaction, 
this phenomenon has been inferred from the expansion of 
a Bose-Einstein condensate in a one-dimensional lattice 
geometry A subsequent experiment Q utilizes the 
reduction of the tunneling matrix element to augment 
the relative importance of interparticle repulsion, such 
that the quantum phase transition from a superfluid to a 
Mott-insulator [H, |j| is induced by adiabatically varying 
the amplitude of the driving force JS- 

These landmark experiments [H, 0] clearly confirm 
that there are efficient control mechanisms for ultra- 
cold atomic gases resulting from time-periodic modu- 
lation. The situation encountered here is akin to the 
dressed-atom approach: An atom in a laser field becomes 
"dressed" by that field and changes its behavior . Sim- 
ilarly, a many-body matter wave becomes dressed in re- 
sponse to time-periodic forcing and acquires properties 
which the unforced, bare matter wave did not have. 

A system of ultracold bosonic atoms in a shaken, suffi- 
ciently deep one-dimensional optical lattice is described, 
in the frame of reference co-moving with the lattice, by 
the driven Bose-Hubbard model defined by the Hamil- 
tonian H(t)^ = H tun + Hint + H dlive (t) @, 0. With b £ 
and fi£ — StS» denoting the bosonic annihilation and 
the number operator for the Wannier state located at 
the site labeled by t = 1,2, ... ,M, one has H tun = 



— J ^fLi 1 {b]bp,+i + ^\-\--\Pe)i wnere the positive hopping 
parameter J implements the kinetics, assumed to be ex- 
hausted by tunneling between adjacent sites. Moreover, 
H mt = f Eti"'("' — 1) with positive interaction pa- 
rameter U describes the repulsion of particles occupying 
the same site. Finally, -ffdrivc(i) = K u cos(wi) J2eLi V-ht 
models time-periodic forcing with amplitude and fre- 
quency to. With the particle number fixed to N, the 
filling n is given by the ratio n = N/M . 

As witnessed by the experiments [H 0], in a time- 



averaged sense the driven system governed by H(t) be- 
haves similar to a system described by the effective, time- 
independent Hamiltonian H Q g = jQ^K^/huj)^^ + ifi„t, 
which means that the effect of the time-periodic force is 
captured by replacing the tunneling matrix element J by 
J c ff = ^(K^/hiv) J, with Jo denoting the ordinary Bessel 
function of order zero. This modification of the hopping 
matrix element is a hallmark of driven quantum tunnel- 
ing [8(; it has been cleanly observed for single-particle 
tunneling in strongly driven double- well potentials 
While it becomes exact for a single particle on a one- 
dimensional lattice endowed with nearest-neighbor cou- 
pling [lflj, the dynamics are considerably more involved 
in the many-body case described by the driven Bose- 
Hubbard model. Due to the manifold ways to create ex- 
citations in the many-body system, the f/eff-description 
is endangered by a multitude of resonances, and holds 
approximately only when Hu is large compared to both 
energy scales which characterize the undriven system, U 
and nJ [H, EH, EH . To further explore the newly emerging 
notion of adiabatic control of driven macroscopic matter 
waves [2j, it is now of great importance to study such 
resonances in detail: When do they occur, how strong 
are they, are they detrimental to coherent control or can 
they, perhaps, even be exploited? These questions mark 
the scope of the present Letter. By means of numeri- 
cal simulations for small systems, we first outline an ex- 
perimentally feasible detection scheme which allows one 
to locate major excitation channels in parameter space, 
and to probe their strengths. We also demonstrate that 
the strength of such excitation channels again is sub- 
ject to coherent control: A resonance can be completely 
quenched by an appropriate choice of the driving ampli- 
tude. In a second step, we make closer contact between 
the dressed-atom picture and the driven matter waves 
considered here by studying their quasienergy spectrum. 
In the final third step we explain our findings quantita- 
tively by means of perturbation theory for Floquet states. 

Consider the following scenario: A system conforming 
to the undriven Bose-Hubbard model H tnn + H mt is pre- 
pared in its ground state for U/J = 0.1. Then a drive 
Hdmvc{t) is switched on, with an amplitude increasing lin- 
early in time, and a high frequency tku/J = 20. Since this 
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FIG. 1: Exact time-evolution of N = 7 particles on M = 7 lat- 
tice sites. Starting in the ground state at interaction strength 
U/J — 0.1, a drive of frequency hcu/J — 20 has been linearly 
ramped up within 50 cycles T — 2tt/ui to the working ampli- 
tude Ku>, before U is increased at various rates n = UT/J = 
0.3 (black), 0.1 (blue), 0.03 (magenta), 0.01 (green), 0.003 
(red), 0.001 (brown). We plot the squared overlap Peff(t) 
of the instantaneous ground state of H e « with the actual 
time-evolved state versus U(t)/ftw at integer t/T. For large 
U(t)/hu, P e ff decreases with decreasing r/. For K^/hu = 2.5 
there is strong resonant excitation at U/hu) — 2/3 (a). For 
K w /huj = 3.4 this resonance is quenched, and another one 
around U/huj — 1 becomes active (b). 



drive is sufficently off-resonant, one expects the system 
to adiabatically follow the ground state of H c g. After the 
working amplitude K u has been reached, it is held con- 
stant. Then the interaction parameter U is ramped up at 
constant rate r\ = UT/J (with T = 2-k/uj) into the regime 
where resonances should make themselves felt. In a lab- 
oratory experiment this can be done, e.g., by increasing 
the transversal confinement used to create the effective 
one-dimensional geometry. We have simulated this pro- 
tocol for a small system with N = M = 7. In Fig. [T] we 
plot the squared overlap P c s{t) = \(ipQ \ip(t))\ 2 of the 
system's true state \tp(t)), obtained by solving the full 
time-dependent Schrodinger equation governed by H(t), 
and the ground states \ipQ^) of the corresponding in- 
stantaneous operators H e s- Figure Q] (a) is obtained for 
K u /fru = 2.5. As expected, P e s stays close to unity 
even when U becomes large, thus validating the H e ff- 
description, until at U/Hlj w 2/3 it decreases suddenly; 
the drop is the more pronounced, the lower the rate r\. 
This abrupt decrease signals resonant excitation. Exper- 
imentally, such resonant excitation can be detected by 



time-of flight absorption imaging. It is indicated by a loss 
of contrast of the sharply peaked structures visible in ei- 
ther the single-particle momentum distribution Q if the 
system is in the superfluid regime (which may be reached 
by a further adiabatic parameter variation), or in the two- 
particle momentum correlations [HI, [3| if the system is 
in the Mott-insulator regime. Interestingly, when choos- 
ing the particular driving amplitude K^/hw = 3.4, this 
excitation channel is closed, and another one at U/hu) ks 1 
appears in Fig. [T] (b). This second resonance is stronger 
than the first one, since the drop is fully developed al- 
ready for larger rj. We conclude: (i) H e a describes the 
system up to surprisingly large interaction strengths U; 
(ii) the excitation observed at U/Hui ss 2/3 in Fig. Q] (a), 
and at U/Hui m 1 in Fig. Q] (b) , cannot be ascribed to a de- 
viation from adiabatic following on the level of H c g, since 
the degree of excitation increases with decreasing param- 
eter variation rate rj; (iii) a resonance can be disabled by 
adjusting the driving amplitude. Thus, by applying this 
or a similar protocol, both the locations and the strengths 
of resonant excitation channels can be probed. 

We now shed light on the physics underlying this de- 
tection scheme, and provide an appropriate theoretical 
framework. Recall that the dressed-atom approach deals 
with atoms interacting with a quantized mode of a ra- 
diation field. Accordingly, the energy level diagram of 
the combined system features identical copies of groups 
of levels displaced against each other by the photon en- 
ergy fku Q. An analogous picture for dressed matter 
waves driven by a classical ti me-p eriodic force is obtained 
by quantum Floquet theory [3, [l6| : Given the Hamil- 
tonian H(t) = H(t + T), one defines the quasienergy 
operator Q = H(t) — ihdt which acts in the product 
space Ji ® T made up from the physical state space 
H and the space T of T-periodic functions, and solves 
the eigenvalue problem Q\u(t))) = e\u(t))). Because 
of the periodic boundary conditions in time, the so- 
lutions have the form \u Vjm (t))) = \u Vj o(t))) exp(imwt), 
with oj = 2tt/T and m = 0, ±1,±2, ...; the label v is 
chosen such that \u„fl(t))) connects to the ^-th energy 
eigenstate when the driving force vanishes. Hence, the 
eigenvalues = e^.o + mhio, called quasienergies, re- 
peat themselves with period hu on the energy axis; each 
state v placing one copy in each "Brillouin zone" of width 
frw. Going back to the actual state space H, the states 
\tp v {t)) = \u Vtm (t))) exp(— i£ v . m t/h) form a complete set 
of solutions to the time-dependent Schrodinger equation. 

Fig. [2] shows a part of the quasienergy spectrum be- 
longing to a small driven Bose-Hubbard system (N = 
M = 5) with fojj/J — 20 and K^/hw — 2 versus U/J. 
Its basic structure, shown in subplot (a), can be under- 
stood as a superposition of copies of the energy spectrum 
of H c ff , shifted against each other by integer multiples of 
huj. The spectrum of H c g possesses bands, made up from 
various types of particle-hole excitations with energies 
roughly corresponding to integer multiples of U, clearly 
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FIG. 2: (a) Quasienergy spectrum of a driven Bose-Hubbard 
system with N — M = 5, hu/J — 20, and K^/huj = 2 ver- 
sus JJ I J. Bands with different slopes belong to different types 
of particle-hole excitations of H c n ■. Resonant coupling of such 
bands results in avoided crossings. The isolated quasienergy 
level, highlighted in (b), emerges from the ground state of the 
undriven system. Clearly visible are the avoided crossings at 
U/fiw « 2/3 (U/J ~ 13) and U/fiw « 1 (U/J « 20) which 
have been detected dynamically in Fig. [T] whereas there are 
no avoided crossings at 1/3 and 1/2. 



identifiable through their slopes. While in Fig. [2] (a) 
quasienergy levels belonging to different copies of the 
^off-spectrum hardly "notice" each other for interaction 
strengths U/J much smaller than hw/ J = 20, there are 
pronounced avoided crossings when U/J becomes compa- 
rable to hw/J, prominently exemplified by the complex 
patterns which appear when U/J is an integer multiple of 
huj/J. Such avoided crossings indicate resonances which 
emerge if eigenstates of H c g are energetically separated 
by an integer multiple of fko] their size quantifies the 
strength of resonant coupling and determines the degree 
of deviation from the i? c ff-description. 

Fig. [2] (b) shows a detail of Fig. [2] (a), focusing on one 
of the quasienergy copies corresponding to the ground 
state of H c ff- After separating from the bands of ex- 
cited states with increasing U/J, thus indicating the 
superfluid-to-Mott insulator transition 0, this level 
crosses several bands associated with different copies of 
the -ffeff-spectrum without being notably affected, until 
it undergoes a wide avoided crossing with such a band 



at U/J Ri ^hu/J w 13, and subsequently an even wider 
one around U/J ~ huj/ J — 20. These avoided crossings 
explain the excitation observed in Fig. Q] The dynamical 
detection scheme illustrated by that figure relies on the 
adiabatic principle for Floquet states [121 ]. With increas- 
ing U, the state \^(t)) adjusts itself to the slowly chang- 
ing parameter and thus follows the quasienergy level cor- 
responding to the ground state of H e g, until it reaches an 
avoided crossing too wide to be passed diabolically. Then 
an incomplete Landau-Zener transition to the anticross- 
ing state excites the system. According to Landau-Zener 
estimates, and in agreemement with the simulations de- 
picted in Fig.[Tl the excitation probability increases expo- 
nentially with both the width of the anticrossing and de- 
creasing parameter speed. Thus, the method of detecting 
resonances in dressed matter waves by monitoring their 
response to slow parameter changes can be regarded as 
a kind of avoided level crossing sprectroscopy. 

Note that in contrast to the regime of linear response, 
suitable for probing properties of the undriven system, 
here we consider the excitation of a system which has 
already been strongly modified by the driving force, in a 
manner described by H c g. Moreover, besides the wide, 
"active" avoided quasienergy crossings there also is a 
host of tiny avoided crossings, reflecting the high den- 
sity of quasienergies in each Brillouin zone, so that ef- 
fectively adiabatic dynamics on the level of H c g actually 
includes fully diabatic Landau-Zener tunneling through 
these narrow anticrossings. In an infinitely large sys- 
tem with a truly dense quasienergy spectrum, the exis- 
tence of a well-defined adiabatic limit cannot, thus, be 
expected 17|. However, realistic parameter variations 
take place on finite time-scales, in all likelihood mak- 
ing the system "blind" against such small features of the 
spectrum. 

We now formalize our reasoning. For each admissi- 
ble set {rif} of site-occupation numbers, we employ the 
usual Fock states |{n^}) = \\i (ntfy~^ (Pi) ne (vacuum) for 
constructing an orthonormal basis of Floquet-Fock states 
\{n e },fh)} = |{n^})exp[-i^sin(wt)^^]exp(imo;t) 
in H <E> T, with fh serving as "photon" index for dis- 
tinguishing different Brillouin zones. Invoking the scalar 
product ((-|-)) = y J^dt (■{■), the quasienergy operator 
Q = Qo + Qi of the driven Bose-Hubbard model pos- 
sesses the matrix elements 

((Wi},m'\Q \{ni},m)} 

= 5fh>,?h({n'(}\(rhtiw + H int + joH tun )\{ne}) , 

((Wt},fh'\Qi\{nt},fh)) 

= (1 - 5fn',7h)js(m-fh')({n'i}\H tVin \{ne}) , 

where j v = ^ v (K^,/hixi) indicates the Bessel function of 
order v evaluated at K^/fuv, and s = X^^f/ 1 ^ — n 'i)^ 
giving s = +1 (s = —1) if fltun tranfers one particle 
by one site to the left (right) [ll|, This splitting of 
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FIG. 3: Coupling strength c 7 of simultaneous resonant exci- 
tation of two particle-hole pairs with the two extra particles 
located at the same site, evaluated at U = 7?lw/3. 

the quasienergy operator is performed such that Qq acts 
within each subspace with fixed "photon" number m in 
a manner conforming to H e g; whereas Q\ describes the 
coupling between these subspaces. 

Let us assume that U is comparable to fuv while hui 
nJ, and treat Qi by perturbation theory. For U 3> n\ J e ff| 
the ground state of H e g is approximately given by the 
extreme Mott-insulator state |MI) = \{ni = n}) with n 
particles localized at each site. Excited states differ from 
| MI) by particle-hole excitations of energy U; these exci- 
tations form bands with widths on the order of ~ n\ J e ff| 
due to tunneling of the particles and holes "on top" of 
| MI). Thus, near U = ahui with integer a = 1,2,... the 
drive is resonant with respect to the creation of a single 
particle-hole pair; eigenstates of Qq differing from |MI) 
by one particle-hole pair and a "photons" are degener- 
ate with |MI) and couple directly (i.e., in first order) via 
Qi by matrix elements of size —yjn(n+ l)Jj s Afh- This 
coupling leads to the large avoided band/level crossings 
visble in Fig. [2] at U/J close to 20 and 40. 

In second order, the simultaneous creation of two 
particle-hole pairs via (quasi-)energetically distant in- 
termediate states is taken into account. Intriguingly, 
second-order coupling between states differing from |MI) 
by (3 "photons" and two separate particle-hole excitations 
of total energy 2U, expected near U = phw/2 with (3 — 
1,3,5 (omitting first-order resonances), vanishes com- 
pletely due to destructive interference between paths in- 
volving different intermediate states. This explains why 
there is no avoided crossing at U/J « 10 in Fig. [2j How- 
ever, there are non- vanishing second-order processes cre- 
ating two overlapping particle-hole pairs, having two par- 
ticles or holes sitting at the same site. Assuming unit fill- 
ing n = 1 , the only possibility is to place both particles at 
the same site, costing the excitation energy 3?7. For such 
excitations near U — jhtu/S with 7 = 1,2,4,5,..., we 
find coupling constants c 7 J 2 n^J (n + l)(n + 2) jfrw with 

strengths c 7 = 5 Eft'=-oo0'(7+m'2^' +3-(i+m')3-m) * 
[(2U/huj — 7 — to') -1 — (U/huj + m') _1 ] which vanish for 
odd 7. The plot of c 7 depicted in Fig. [3] testifies that 
these strengths depend in an oscillating manner on the 
driving amplitude. In particular, it is possible to adjust 



that amplitude such that the resonant coupling strength 
vanishes. For instance, the zero of ci at K^/hw w 3.4 
is the reason for the resonance quenching illustrated in 
Fig. [TJ In vtYi order, coupling matrix elements generally 
are ~ nJ(nJ/hui)"~ 1 ; however, we have hardly noticed 
third-order effects in our numerical simulations. Thus, 
degenerate-state perturbation theory in H (g> T system- 
atically uncovers the hierarchy of resonances which, in a 
system with slowly changing parameters, become observ- 
able order by order with decreasing parameter speed. 

To conclude, we have outlined a scheme for probing res- 
onances which endanger the adiabatic control of macro- 
scopic matter waves achievable through time-periodic 
forcing 0]. The theoretical analysis of this scheme re- 
veals far-reaching conceptual similarities between dressed 
atoms and dressed matter waves in shaken optical lat- 
tices, thus opening up wide new grounds between quan- 
tum optics and matter- wave physics. 
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